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TWO-DIMENSIONAL DIFFUSION THEORY ANALYSIS OF REACTIVITY 

EFFECTS OF A FUEL-PLATE-REMOVAL EXPERIMENT 

By Edward R. Gotsky^ James P. Cusick, 
and Donald Bogart 


SUMMARY 

Two-dimensional two-group diffusion calculations were performed on 
the NACA reactor simulator in order to evaluate the reactivity effects 
of fuel plates removed successively from the center experimental fuel 
element of a seven- hy three-element core loading at the Oak Ridge Bulk 
Shielding Facility. The reactivity calculations were performed hy two 
methods: In the first, the slowing-down properties of the experimental 

fuel element were represented hy infinite media parameters; and, in the 
second, the finite size of the experimental fuel element was recognized, 
and the slowing-down properties of the surrounding core were attributed 
to this small region. The two calculation methods agreed reasonably 
well with the experimental reactivity effects. 


INTRODUCTION 

The NACA Lewis laboratory is interested in a high-flux research re- 
actor with which the components of high power density reactors under de- 
velopment for various flight applications can be studied. The reactor 
being considered employs an array of aluminum fuel elements similar to 
those of the Materials Testing Reactor in Idaho. 

The geometric center of the active lattice of a reactor is an at- 
tractive location for experiments because of high flux and flux symmetry. 
The geometric center of the core, however, is also the most sensitive 
reactivity region for accidental compositional cheinges, and great care 
must be exercised in the design of such in-pile experiments. Many re- 
activity calculations have been made for both fueled and unfueled exper- 
iments in a center test hole in support of the design of the NACA re- 
search reactor. Since reactivity calculations are subject to uncertainty, 
an experimental program with the Bulk Shielding Reactor at Oak Ridge 
National Laboratory was initiated in order to evaluate the methods of 
analysis. In these Oak Ridge experiments, the core configuration of the 
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MCA research reactor was mocked up within the limits of excess reactiv- 
ity and materials available to the Bulk Shielding Reactor. The loading 
was such that the reactivity effects of voids within the core could be 
measured and the worth of fuel-element water passages and fuel-element 
plates could be determined. The Bulk Shielding Reactor is fully dis- 
cussed in reference 1. The mipublished experimental data for this par- 
ticular loading were obtained at the Oak Ridge National Laboratory by 
E. B. Johnson, K. M. Heniy, J. D. Kingdon, and T. M. Hallman. 

The experiments in which the reactivity effects are measured are 
conveniently analyzed by group diffusion calculations. Within the limits 
of diffusion theory, two- and three-group analyses yield reasonable neu- 
tron flux distributions for thermal and epithermal reactors. For 
hydrogen-moderated reactors, many modified group diffusion methods have 
been used to interpret criticality experiments (refs. 2, 3, and 4). 

The present paper is concerned with the Bulk Shielding Reactor fuel- 
plate-removal experiments and the correlation of these data with diffusion 
theory calciilations . These experiments have been analyzed by two-group 
two-dimensional diffusion calculations. The solutions have been obtained 
on a two-dimensional nuclear-reactor simulator that is based on the de- 
sign of the simulator described in reference 5. 

The cooperation of the Oak Ridge National Laboratory and the staff 
of the Bulk Shielding Facility in performing the reactivity experiments 
is gratefully acknowledged. The authors are greatly indebted to Miss 
Dorothy Hood, who supervised the calculations on the reactor simulator. 


RMCTOR AND REACTIVITY EXPERIMENTS 

The Bulk Shielding Reactor is an assembly of fuel elements that may 
be arranged into various critical configurations. The fuel is highly 
enriched uranium contained in aluminum-clad fuel plates. A complete fuel 
element is made up of 18 fuel plates and contains a total of about 140 
grams of uranium- 235. 


Reactor Loading 

The reactor configuration used in the reactivity experiments is 
shown schematically in figure 1. The loading consisted of 21 fuel ele- 
ments in a seven-by-three array with beryllium oxide (BeO) reflector 
pieces arranged in one row on the north and two rows on the south. The 
east and west faces of the core were reflected by water and permitted 
complete insertion or complete withdrawal of the two guillotine safety 
blades® Each safety blade consisted of a thin cadmium sheet between two 
aluminum plates approximately 12 by 24 inches . The guillotine safety 


099b 


CJ-1 back 


NACA TN 4164 


3 


b^des were guided and positioned^ relative to the core, by grooved alu- 
minum pieces and were supported by electromagnets actuated by scram cir- 
^ cults. A regulating rod and a core safety rod were provided in the 

control-rod fuel elements in grid positions 23 and 27, respectively. 

The regulating-rod calibration was obtained by the method of dis- 
tributed poisons for three configurations of the guillotine safety 
blades; namely, both guillotines out, number 1 guillotine out with 
number 2 guillotine in, and both guillotines in. The elements through 
o conventional safety rod and the regulating rod move contained 

half the normal number of fuel plates and, therefore, about 70 grams of 
^anium-235 each. The primary-reflector elements were hot-pressed 
eryllii^ oxide blocks encased in watertight aluminum cans of the same 
outer dimensions as the fuel elements. The reactor was moderated and 
cooled by water that also served as the secondary reflector and reactor 
shield. center grid position 25 was occupied by the experimental 


It should be noted that this core loading had vertical symmetry 
about the horizontal midplane and, in addition, east-west symmetry about 
a vertical plane across the narrow core dimension. (Although the top 
and bottom reflectors were not identical, the reactor was symmetrical 
nuclearwise because the top and bottom reflectors were predominantly 
^ter.) The vertical symmetry permits a typical horizontal slice through 
the reactor to be analyzed. Because of symmetiy within the plane of 
t^s horizontal slice, only one-half of this slice must be simulated for 
diffusion calciilations . An additional calculation indicated that the 

partially inserted core control rod does not significantly 
affect the flux distributions at the center grid position (25). 


Removable-Fuel-Plate Assembly 

In the reactivity experiments discussed in this report, grid posi- 
tion 25 was occupied by a standard fuel element in which only 6 of the 
18 c^ed fuel plates were brazed into the grooved side plates. A photo- 
graph of this removable-fuel-plate assembly is shown in figure 2. The 
fuel plates are removable with the exception of the two end plates and 
the four plates blocked by the fuel-element handle. The plates are 60 
mils thick, and the water passages are 117 mils thickj the fuel element 
IS 24 inches long and 3 inches wide. 

These fuel plates were removed in two sequences. In sequence I 
the interactions were maximized by successively removing adjacent fuel 
p^tes and thereby foiming an enlarging water region near the center of 
the core. In sequence II, the interactions of each plate removed were 
initially minimized by removing plates farthest apart from each other. 
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These seq,uences are illustrated in figiire Z, in which letters A to L 
represent the fuel plates removed in each step. 


Experimental Results 

The reactivity per incremental fuel plate removed in the maximum 
interaction (sequence l) and in the minimum interaction (sequence II) is 
presented in figure 4. The net effect of the removal of fuel plates and 
the consequent enlargement of water regions decreased the reactivity in 
all cases. The reactivity for each fuel plate removed was greatly af- 
fected by the sequence of removal. In sequence I, for example, the 
worth of adjacent fuel plates increased rapidly as the region devoid of 
fuel enlarged. On the other hand, in sequence II, individual fuel plates 
from opposite sides of the element had equal worth, while the worth of 
fuel plates from intermediate positions in the element increased slowly 
at first. 

In these removable-fuel-plate sequences, available core excess re- 
activity permitted 11 fuel plates to be removed. In the maximum inter- 
action, removal of the twelfth plate (number D in fig. 4(a)) shut down 
the reactor when the plate was 45-percent withdrawn from the core. At 
this point, the regulating rod was fully withdrawn from the core; this 
indicated that the withdrawn portion of the twelfth fuel plate was worth 
0.16 AK/K percent. (Symbols are defined in appendix A.) The accumu- 
lated reactivity of the first 11 fuel plates was 2.36 percent, 'i^e com- 
plete twelfth fuel plate is estimated to be worth at least 0.16/0.45 = 
0.36 AK/k percent (shown in fig. 4(a) as a tailed data point). In 
the minimum interaction, the worth of the twelfth fuel plate is readily 
estimated by subtracting the cumulative reactivity of sequence II, step 
11, from the cumulative reactivity of sequence I, step 12. A reactivity 
of 0.48 Z^/k percent was obtained, as shown in figure 4(b) as a tailed 
data point. 

The cumulative reactivity of these sequences is presented in fig- 
ure 5 as a function of the number of fuel plates removed (each fuel 
plate contains 7.78 g of uranium- 235 ) . The minimum-interaction sequence 
was used as a basis for the reactivity calculations reported herein and 
for the comparison between the analytical group diffusion results and 
the experimental data. 

GROUP DIFFUSION CALCULATIONS ON REACTOR SIMULATOR 

The Simulator 

The MCA two-dimensional reactor simulator is an analogue computer 
for obtaining solutions to the neutron-group diffusion eqmtions. 
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Reactivity problems in which spatial symmetry exists in only one dimen- 
sion may be solved. The simulator is a direct extension of the one- 
dimensional device in use at the NACA Lewis laboratory since 1951 
(ref. 5). 

The two-dimensional reactor simulator is shown in figure 6. The 
fast- and thermal- neutron- group resistor-network boards^ in which the 
space-point resistors for a particular problem are prefabricated in 
aluminum plug-in units^ may be seen on the right. The source-input 
panels are shown on the left. The problem is iterated from group to 
group by varying the source input until a converged flux distribution 
and value of the effective multiplication factor Keff are obtained. 

In the two-group formulation used^ epithermal fission and absorption are 
part of the fast-group processes. 


Simulator Representation of Loading 

The core loading shown previously had vertical symmetry about the 
horizontal midplane and^ in addition^ symmetry about a vertical plane 
across the narrow core dimension. The vertical symmetry permits a typi- 
cal horizontal slice through the reactor to be analyzed. Because of 
symmetry within the plane of this horizontal slice^ it is necessaiy to 
simulate only half of this slice for the diffusion calculations. Repre- 
sentation of this geometry on the reactor simulator^ therefore^ takes 
the fonn shown in figure 1 y where the core fuel elements^ the primary 
beryllium oxide reflectors^ and the secondary water reflector are viewed 
from above. The fuel element in the central core position^ from which 
fuel plates were successively removed^ is shown in the inset. There is 
a net point for each area shown. Each coarse net point has the dimen- 
sions of a fuel element^ about 3 inches square. The medium- sized net 
points are one-fourth of the coarse grid size. The area at the center 
of the core represents one-half of the experimental fuel element from 
which fuel plates were individually removed. In order to obtain a de- 
tailed description of the fluxes in this region^ fine-sized net points 
that represent areas about a centimeter on a side were used. In all^ 
the calculation involved 69 coarse net points^ 34 medium net points, and 
32 fine net points in each group of the two-group solution. The nuclear 
properties used for each space point are for the homogenized contents of 
the experimental region represented. 


Two -Group Formulation 

The two-group multiregion flux equations and a brief description 
of the methods of evaluating group parameters are given in appendix B. 
The two-group constants used in these calculations for the unperturbed 
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core are shown in table I. These parameter values are similar to those 
used to predict reactivity effects in the early Materials Testing Reactor 
(Idaho) loadings, where fuel elements similar to those in the Bulk 
Shielding Reactor were used. In the two-group formulation used here, 
epithermal fission and absorption are retained as part of the fast-group 
processes; this requires the use of the fast multiplication constant 

and the nonabsorption probability P-th- shown in reference 6 

that a two-group formulation including epithermal absorption and fission 
as used herein predicts reactivities that agree with a formulation which 
considers in detail the energy-dependent nuclear processes in the epi- 

2 2 

thermal region. Values of and derived from experimental 

values for metal-water mixtures; and other nonthermal values are calculated 
from slowing-down distributions in infinite media regions of the same com- 
position. Parameter values. for regions in the center fuel element with 
fuel plates removed were estimated in this manner also. 

In all cases, the neutron leakage in the direction nonnal to the 
plane of the two-dimensional solution is assumed proportional to the 
local flux with the vertical geometric buckling as the constant of 

proportionality. The b| value of 0,001622 is based on an actual core 

height of 62.0 centimeters plus reflector savings of 16.6 centimeters 
due to the aluminum- water regions above and below the core. The problem 
of keeping this vertical buckling constant as more and more fuel plates 
are removed accounts for part of the small discrepancy between exper- 
iment and calculations. 

The question immediately arises as to the validity of diffusion 
theory in regions that are geometrically small, such as regions where a 
few fuel plates have been removed. For example, removal of two adjacent 
fuel plates results in a water passage that is about 1 by 7 centimeters 
in cross section and runs the full core height. The slowing-down proper- 
ties of such a small region can hardly be characterized by the slowing- 
down length in water, which is about 5.6 centimeters. It is reasonable 
to assume that the slowing-down properties of the surrounding core may 
be more significant. For this reason, the diffusion calculations for 
each core configuration were perfoimed twice; first, using the homogeneous 
diffusion properties for each region regardless of size (referred to 
hereafter as ’’local region slowing-down” cases); and second, using the 
values of and p^^ for the surrounding core to characterize the 

slowing down with absorption in the regions where fuel plates were re- 
moved (referred to as ’’surrounding region slowing-down” cases). In the 
surrounding region slowing-down cases, two-group constants other than 
and p^j^ were based on the local region properties. 
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Experimental Configurations Calculated 

Calculations of the reactivity effects resulting from fuel-plate 
removal were performed for the reference core with no fuel plates re- 
moved and for the minimum interaction sequence in which 6, 10, and the 
equivalent of 13.5 fuel plates were removed from the center fuel element. 
The fine-grid-point arrangement used in each case duplicated composi- 
tionally the experimental geometry within the limits of the 32 space 
points available. An additional case, in which the range of the experi- 
mental data was extended by removing all 18 plates of the central fuel 
element, was calculated in order to estimate the worth of one complete 
fuel element at the center of the core. 


CALCULATED RESULTS 
Unperturbed Core 

The representative two-dimensional flux distributions are shown in 
figure 8 for the unperturbed core with no fuel plates removed from the 
center fuel element. Figure 8(a) shows the integrated fast-flux distri- 
bution in the core and reflectors. The magnitude of the fluxes is shown 
relative to a spatial average core thermal-neutron flux of unity. Fig- 
ure 8(b) shows the thermal-neutron flux distribution with characteristic 
peaking in the reflectors and falloff toward the corners of the slab. 
There is an over-all asymmetry due to the thicker primary reflector on 
one face of the core. 

To observe the effects of fuel-plate removal from the central fuel 
element, the flux distributions across the narrow core dimension where 
maximum perturbations occur are noted. Reactivity effects were evaluated 
by integrating the two-dimensional flux distributions obtained over the 
core volume. 

On an absolute basis, the value of Keff for this Bulk Shielding 
Reactor core loading was observed experimentally to be 1.025; that is, 
the loading had 2.5-percent excess reactivity. The value of Keff ’ 

estimated from two-dimensional diffusion calculations is 1.032, an esti- 
mated excess reactivity of 3.2 percent. 


Local Region Slowing-Down Model 

Neutron flux distributions for several local region slowlng-down 
cases are shown in figure 9. The flux distributions across the narrow 
core dimension are shown for three cases: (l) the unperturbed core, 

(2) 10 fuel plates removed, and (3) all 18 fuel plates in center fuel 
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element removed. The fluxes are relative to a spatial average thermal 
flux of unity in the core. Shown in figure 9 are data for the core with 
the central fuel element^ the beryllium oxide side reflectors of differ- 
ent thicknesses^ and the secondary water reflectors. Removal of the fuel 
plates and^ therefore^ fast-neutron sources^ with consequent replacement 
by water^ has reduced the fast flux (p ^ for a large region in the vicin- 
ity of the central fuel element by providing a large internal slowing- 
down region. The thermal fluxes cpth consequently built up in this 

region. The overflow of neutrons diffusing into surrounding fuel ele- 
ments is felt for about two thermal- diffusion lengths into the core, 
after which the thermal fluxes are unperturbed. The reactivity worth of 
remaining center fuel plates and of plates in adjoining elements is in- 
creased as a result. The net reactivity effect is the sum of the -AK 
due to fuel-plate removal and the -iAK due to adjoining fuel plates op- 
erating in higher thermal fliixes. 


Surrounding Region Slowing-Down Model 

Figure 10 illustrates the two-group neutron fluxes for 10 fuel 
plates removed^ as calculated by assuming the values of L^* and p^^^ 

for the missing-fuel-plate region to be the same as those of the sur- 
rounding core. These flux distributions are compared with the fluxes for 
the unperturbed core. The slowing-down parameter for the core is 

p 1 

larger than the value of Lf for the missing-fuel-plate region, with the 

result that fewer neutrons slow down in this central region. Consequent- 
ly, the thermal-neutron flux at the center region for this surrounding 
region slowing-down case is lower than the thermal-neutron flux for the 
corresponding local region slowing-down case. The remaining fuel plates 
in the center test element and in the adjoining elements, being in re- 
gions of lower thermal flux for the surrounding region slowing-down case, 
are worth less than the fuel plates for the local region slowing-down 
case, where thermal flux is higher. Inasmuch as the net reactivity ef- 
fect is the Slim of the -AK due to fuel-plate removal and the -h^K due 
to adjoining fuel plates operating in higher relative thermal fluxes, 
quite different effects may be expected from the two diffusion 
calculations. 


DISCUSSION 

The net reactivity effects of fuel-plate removal from the center 
fuel element, as measured at the Bulk. Shielding Reactor and as presently 
calculated, are compared in figure 11. The local region slowing-down cal- 
culations consistently underestimate the observed reactivity effects. On 
the other hand, the surroxmding region slowing-down calculations predict 
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larger reactivity effects for fuel-plate removal up to about 10 plates, 
so that the actual physical situation is bracketed. 

Both calculations underestimate an extrapolation of the experi- 
mental data corresponding to removal of the remaining fuel plates in 
the center fuel element. This may indicate that a third-dimensional ef- 
fect is becoming important. In all of the present calciiLations, neutron 
leakage in the vertical direction, tliat is, normal to the plane of solu- 
tion, was taken to be proportional to the local flxnc <P with a fixed 
value of the vertical geometric buckling as the constant of propor- 

tionality. This assumes a cosine flux distribution in the z-direction 
above and below the plane of solution. Furthermore, this assumption im- 
plies that the flux in any horizontal plane is directly proportional to 
the flux in the horizontal midplane. 

As the mis sing- fuel-plate region at the center of the core enlarges, 
however, an increasingly important diffusion region with reflector proper- 
ties is provided^ this gives rise to a nonuniform flux distribution above 
and below the plane of solution. This internal water region acts essen- 
tially as a fast control rod by collecting fast neutrons over a relative- 
ly l^'i’ge range of several fast diffusion lengths in the core, by moder- 
atlng, and then by diffusing them back into the core as thermal neutrons 
over a shorter range of several thermal diffusion lengths. 

To illustrate this third-dimensional effect, consider figure 12, in 
which portions of the center fuel element are removed from the top and 
bottom of the core and replaced with water. This internal reflector 
(water) gathers in fast neutrons from all directions over a range of 
about one fast-transport mean free path in the core (4 cm) and moves the 
neutrons away from the center of the core. These neutrons are therma- 
lized in the water region and diffuse into the core over a range of about 
one thermal-transport mean free path (0.8 cm). The net result is a trans- 
fer of fast neutrons away from the center of the core to regions of less 
iniportance. As the internal-reflector region increases in height, the 
migration of fast neutrons away from the core center also increases. 

The two-dimensional diffusion calculations cannot account for this 
fast-neutron migration away from the center of the core, and therefore 
higher thermal fluxes are indicated in the central core region than are 
present in the experiment. The reactivity worth of the fuel elements 
surrounding the central fuel element are consequently overestimated, and 
the worth of the central fuel element is underestimated. 


An approximation of the diffusion in the internal water region can 
be made by decreasing the effective average height of the core; this re- 
sults in an increase in and a consequent increase in the vertical 

leakage term DfBz. To illustrate this, a reduction in effective core 
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height of 1 centimeter results in a decrease in reactivity of 0.3 percent 
/iK for this core loading. (The reduction in reactivity is linear for 
about 6 cm in effective core height.) 

CONCLUSION 

The three-dimensional reactivity effects measured in the Bulk Shield- 
ing Reactor core-loading experiments are approximated by a two-dimensional 
calculation. Over the range of experimental data^ the two-dimensional 
diffusion calculations predict reactivity effects that bracket the ex- 
perimental data reasonably well. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland^ Ohio^ October 18, 1957 
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APPENDIX A 
SYMBOLS 



buckling in direction normal to plane of solution 

D 

o 

CC 

^ h 

diffusion coefficient 
height of active lattice 

K 

multiplication constant 

Keff 

effective multiplication factor for reactor 


mean-square slowing-down length for fission neutrons 
mean- square diffusion length for thermal neutrons 

1 ■. 

fraction of absorbed neutrons escaping capture in slowing down 

D 

^ 

top reflector savings 

^2 

bottom reflector savings 


macroscopic absorption cross section 


macroscopic fission cross section 

^S1 

macroscopic slowing-down cross section 

V 

nimiber of neutrons emitted per fission 

(P 

neutron flux 


Subscripts : 


f 

fast neutron group 

th 

thermal neutron group 


12 


MCA 03J 4164 


APPENDIX B 


DIFFUSION EQUATIONS AND GROUP CONSTANTS USED IN 
REACTIVITY CALCULATIONS 
Equations 

Two-group multiregion two-dimensional diffusion calculations are 
based on the following flux equations that include effects of fast ab- 
sorption and fission: 


/ \ 

/ ^ 2 ^2 \ 

®th ( |j2 + ‘‘^th - (^A,th + D^-hB2)^th + = 0 (B2) 


In equations (B1) and (B2), effective absorbers are introduced to ac- 
count for neutron leakage in the direction normal to the plane of solu- 
tion; these are D^B| and the epithermal and thermal groups, 

respectively. The quantity B^ is the geometric buckling given in terms 

of the equivalent bare-core dimensions in the direction normal to the 
plane of solution. For the present reactor. 


B 


2 

z 


(h + + 62)' 


Inasmuch as the macroscopic cross sections are the effective averages 

in each energy group, the mean- square slowing- down length L§ and the 

2 ^ 
mean-sqiiare thermal diffusion length Lth may be defined as 


i|=== 

^A,f ^Sl,f 


(B3) 


t2 - 

^th - ? 


■"th 


A,th 


(B4) 
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The fraction of absorbed neutrons escaping capture in slowing down 
^•Iso may be expressed in terms of these effective values as 


2si,f 

^A,f %l,f 

(B5) 

The multiplication constants representing the 
born per neutron absorbed in each group are 

number of neutrons 

^A,f 

(B6) 

^ %’,th 

(B7) 


With these definitions, two-group equations (Bl) and (B2) maybe 
rewritten as : 



®th 



+ 



‘Pth - 


Pth 

4h 






Pth ,2 ‘Pf 


Df 

^ (D^ = 0 


(B8) 

(B9) 


Equations (B8) and (B9) will apply for passive regions if 
% = =0. At the interfaces between regions, the fluxes and neutron 

currents are continuous. At the extrapolated boundaries of the secondary 
reflector^ the fluxes are zero- 


Group Parameters 

The constants required for solution of the two-group equations (B8) 
and (B9) in each region are: 

(1) Fast: l2, d^, K^, p^^ 

(2) Iliermal: 

The validity of a group representation depends greatly on the cor- 
rectness of the group constants used to represent the various competing 
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nuclear processes . Procedures used to evaluate group constants for 
highly enriched uranium hydrogeneous media are described in references 4 
and 6 and are outlined in the following section. 

Fast constants . - For neutrons slowing down in hydrogeneous media, 
the bulk of the contribution to is made at high neutron energies, 

that is, at energies above the epithermal region for which the inverse 
velocity variation of microscopic absorption cross section has reduced 
absorption to a negligible value. For highly enriched uranium systems, 
most of the absorption occurs in the epithermal region in the energy 
range below about 1000 electron volts. Hence, fast-group diffusion, 
leakage, and slowing-down processes may be effectively separated from 
epithermal absorption processes. If Fermi age theory is taken to apply 
in the epithermal region, the residual slowing-down, diffusion, and 
absorption processes in this region may be treated in detail. To ac- 
complish this, neutrons should be divided into three separate groups: 
fast, epithermal, and thermal, as in References 4 and 6. All neutron ab- 
sorption, and therefore production, is restricted to the epithermal and 
thermal regions . Group-diffusion slowing down in the fast region is then 

characterized by values of l| estimated from correlations of the few ex- 
perimental metal-water-mixture measurements (ref. 4). Slowing down in 
the epithermal region is calculated by Fermi age theory, which provides 
the variation of neutron flux with energy to permit evaluation of the 
effective values of 2^, 2p, and p. These average values are obtained 

by weighting local values by the neutron flux in the epithermal region. 

In the present two-group calculations, epithermal fissions and 
absorptions have been retained in the fast group, and the age-theory 
neutron slowing-down distribution in infinite media has been used to 
evaluate the fast-group constants. The fast constants D^, Kf, and p^j^ 

are, therefore, obtained by weighting local values according to the 
energy distribution of neutron flux in an infinite medium of the same 
composition as the pertinent reactor region and as indicated by age 
theory. The fission spectrum is included in determining the energy 
distribution of neutron flux. 

It is shown in reference 6 that a two-group formulation including 
epithermal absorption and fission predicts reactivities that agree with a 
formulation which considers in detail the energy -dependent nuclear pro- 
cesses in the epithermal region. 

Thermal constants . - The thermal parameters are calculated from 
equations (B4) and (B7) . Values of 2^ and ,th macroscopic 

cross sections obtained by averaging local values over the Maxwellian 
distribution of thermal neutron flux. The effects of thermal spectrum 
hardening due to preferential absorption in the lower energy portions of 
the thermal distribution have been neglected. 
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The diffusion coefficient is given by the sum of the contri- 

butions of each constituent atom, the hydrogen atom usually being the 
largest contributor. The effects of the chemical binding of the hydro- 
gen atom to the water molecule, however, significantly alter thermal 
scattering properties of hydrogen because of the increasing effective 
mass of the hydrogen atom for neutrons in the lower energy portions of 
the thermal distribution. 

The effective value of Dth for any media is also obtained by av- 
eraging local values over the Maxwellian distribution of neutron flux, 
provided the variation of effective mass of the hydrogen atom with lo- 
cal energy is known. In the absence of these specific data, a method 
attributable to Radkowsky (described in ref. 7), which checks experi- 
mentally determined thermal diffusion properties of water for a range 
of temperature, was used. A value of 31.4 barns was obtained for the 
effective microscopic transport cross section for hydrogen. This value, 
which corresponds to an experimental value of thermal diffusion length 
of 2,80 centimeters in water at room temperature, has been used in the 
present calculations. 
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TABLE I. - TWO-GROUP DIFFUSION PARAMETERS 


[Vertical buckling in all cases, = 0.001622.] 


Eeactor region 

Fast group 

Thermal group 

Df 

l| 

Kf 

Pth 

®th 

Lth 

^th 

Bulk Shielding Reactor 

1.320 

64.0 

1.517 

0.864 

0.269 

3.821 

1.612 

core (0.423 Al, 








0.577 H 2 O, U^^^) 








BeO reflector 

.783 

117.5 


1.000 

.563 

242.0 


(0.876 BeO, .074 Al, 








.050 H 2 O) 








HgO reflector 

1.146 

31.4 


.977 

.160 

8.189 




Beryllium oxide reflector 


Standard fuel assembly 


North 


Experimental fuel assembly 


Reactor grid plate 


Number 2 guillotine 
safety blade 


' — Number 1 safety rod 


Regulating rod 


Number 1 guillotine safety blade 


Figure 1. - Reactor loading for reactivity experiments at Bulk Shielding Reactor. 
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Figure 2. - EemoTable-fuel-plate element. 
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Standard fuel element 



12 Removable fuel plates, ABODE FGHIJKL 


Maximum Interaction (sequence l) 
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Figure 3. - Sequences of removing fuel plates In standard fuel element of 
Bulk Shielding Reactor. 
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Figure 4. - Reactivity for each fuel plate removed. 
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Figure 5. - Cumulative reactivity of fuel plates in element at grid position 25. 
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Figure 6. - Two-dimensional reactor simulator 
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Removable -fuel-plate assembly at center of core 


Figure 7. - Reactor geometry for two-dimensional simulator solutions. 
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(a) Fast flux. 

Figure 8. - Flux distribution for unperturbed core. 
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Reactor length, in. 
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(h) Thermal flux. 

Figure 8. - Concluded. Fliax distribution for unperturbed core. 
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Figure 9 . - Fl\ix distributions along reactor width for local si, wing-down cases 
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Figure 10. - Flux distributions along reactor width; slowing-down 
comparisons for cases with 10 plates removed. 


MCA TN 4164 



Figure 11. - Comparison of experimental and calculated reactivity effects of fuel-plate 
removal . 
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Figure 12. - Illustration of asymmetry due to water passage in core center. 
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